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ABSTRACT 

We measure the surface density of z'-band dropout galaxies at z ~ 6 through wide-field HST/ ACS 
imaging and ultra-deep Keck/DEIMOS spectroscopy. Using deep EST/ ACS SDSS-i' (F775W) and 
SDSS-z' (F850LP) imaging from the Great Observatories Origins Deep Survey North (GOODS-N; 
200arcmin 2 ), we identify 9 i'-drops satisfying {%' — z')ab > 1-5 to a depth of z' AB — 25.6 (correspond- 
ing to T\jy at z ~ 3). We use HK' imaging data to improve the fidelity of our sample, discriminating 
against lower redshift red galaxies and cool Galactic stars. Three i'-drops are consistent with M/L/T 
dwarf stars. We present ultra-deep Keck/DEIMOS spectroscopy of 10 objects from our combined 
GOODS-N and GOODS-S i'-drop sample. We detect Lyman-a emission at z = 5.83 from one object 
in the GOODS-S field, which lies only 8arcmin (i.e. 3 h^} Mpc) away from a previously confirmed 
z = 5.78 object. One possible Lyman-a emitter at z = 6.24 is found in the GOODS-N field (although 
identification of this spatially-offset emission line is ambiguous). Using the rest-frame UV continuum 
from our 6 candidate z ~ 6 galaxies from the GOODS-N field, we determine a lower limit to the unob- 
scured volume-averaged global star formation rate at z ~ 6 of (5.4 ± 2.2) x 10~ 4 ft.70 M Q yr _1 Mpc -3 . 
We find that the cosmic star formation density in Lyman Break galaxies (LBGs) with unobscured star 
formation rates > 15-Moyr" 1 falls by a factor of 8 between z ~ 3 and z ~ 6. Hence the luminosity 
function of LBGs must evolve in this redshift interval: a constant integrated star formation density at 
z > 3 requires a much steeper faint-end slope, or a brighter characteristic luminosity. This result is in 
agreement with our previous measurement from the GOODS-S field, indicating that cosmic variance 
is not a dominant source of uncertainty. 

Subject headings: galaxies: formation - galaxies: evolution - galaxies: starburst - galaxies: high 
redshift - ultraviolet: galaxies - surveys 



1. INTRODUCTION 

In recent years the identification and study of very high 
redshift galaxies (z > 5) has been an area of active and 
rapidly-advancing research. The increasing availability 
of 8- and 10- meter class telescopes, combined with the 
development of modern instrumentation such as the Ad- 
vanced Camera for Surveys (ACS) on the Hubble Space 
Telescope (HST) has allowed imaging of ever fainter and 
more distant galaxies. 

While large ground-based telescopes have made pos- 
sible spectroscopic confirmation of rcdshifts for galax- 
ies beyond z ~ 6 (e.g. Hu et al. 2002), this process 
is expensive in telescope time and is only possible for 
those objects with strong emission lines or which are 
lensed by intervening objects (e.g. Ellis et al. 2001). 
As a result, the use of "photometric redshift" selection 
from broad-band colors (such as the Lyman-break tech- 
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nique of Steidel, Pettini & Hamilton 1995) has become a 
widespread method of identifying large samples of high- 
redshift galaxies. 

In Stanway, Bunker and McMahon (2003, hereafter pa- 
per 1) we extended the Lyman-break selection technique 
to z ~ 6, and described the photometric selection of 
i'-band drop-out candidate high redshift objects in the 
Chandra Deep Field South (CDFS). We demonstrated 
that public data from the Advanced Camera for Surveys 
(ACS - Ford et al. 2002) on HST, released as part of 
the Great Observatories Origins Deep Survey (GOODS - 
Dickinson & Giavalisco 2003 1 ), is of sufficient depth and 
volume to detect very high redshift galaxies, and that 
color selection can be used to reject less distant objects. 
We used this data set to examine the space density of 
UV-luminous starburst galaxies at z w 6 in the GOODS- 
South field (the CDFS) and to place a lower limit on the 
comoving global star formation rate at this epoch. Fur- 
thermore, in Bunker et al. (2003, hereafter paper 2) we 

1 see http://www.stsci.edu/ftp/science/goods/ 
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presented deep spectroscopy for object SBM03#3, con- 
firming its redshift as z — 5.78 and illustrating the effec- 
tiveness of our i'-drop selection. 

In Section[5]of this paper we present a similar selection 
of 9 i'-drop objects with large ii! — z')ab colors in the 
GOODS/4 OS northern field (hereafter GOODS-N), cen- 
tered on the Hubble Deep Field North (HDFN, Williams 
et al. 1996). Deep spectroscopy using the DEIMOS spec- 
trograph on the Keck II telescope is presented in Sec- 
tion[3]for half of our (i 1 — z')ab > 1-5 sample, including 
5 i'-drops in the GOODS-N. We also revisit our GOODS- 
S i'-drop sample from paper 1, and present 4 spectra in 
this field (in addition to the z = 5.78 GOODS-S galaxy 
reported in Paper 2). 

In Section^] we discuss the effects of completeness cor- 
rections and surface brightness effects on our sample. In 
Section [3] we estimate the global star formation rate at 
z ~6, and contrast the GOODS-N dataset with that for 
the GOODS-S, enabling us to address the issue of cos- 
mic variance. We also discuss the effects of completeness 
corrections, dust extinction and evolution of the Lyman- 
break galaxy luminosity function, and we compare our 
results to other recent work in this field. Our conclu- 
sions are presented in Section HO 

For consistency with our earlier work, we adopt the 
following cosmology: a A-dominated, flat Universe with 
0\ = 0.7, ft M = 0.3 and H = 70 /i 7 okms _1 Mpc" 1 . All 
magnitudes in this paper are quoted in the AB system 
(Oke & Gunn 1983) and the Madau (1995) prescription, 
extended to z = 6, is used where necessary to estimate 
absorption due to the intergalactic medium. 

2. GOODS-NORTH /-DROPS: PHOTOMETRY AND 
CANDIDATE SELECTION 

2.1. GOODS HST/ACS imaging and i' -Drop Selection 

We present here a high redshift candidate selection 
based on the first four epochs of GOODS HST/ACS ob- 
servations of the GOODS-N available at the time of writ- 
ing. Each epoch of data comprises imaging of either 15 or 
16 adjacent 'tiles' in each of the F606W (V, 0.5 orbits), 
F775W (SDSS-i', 0.5 orbits) and F850LP (SDSS-z', 1.0 
orbits) broad-band filters. Three orbits in the F435W 
i?-filter were obtained in a single epoch at the start of 
the observing campaign. As in paper 1, we analyze v0.5 
of the publically-released reduced data 2 . 

We use the photometric calibration determined by the 
GOODS team: 

ma §AB — zeropoint — 2.5 log 10 (Count rate/s _1 ), 

and note that the zeropoints for the ACS instrument 
changed between epochs 1 and 2 of the GOODS-N obser- 
vations. In epoch 1 AB magnitude zeropoints for the v- 
band (F606W), i'-band (F775W) and z'-band (F850LP) 
are 26.505, 25.656 and 24.916 respectively. In Epoch 2 
and thereafter the zeropoints for the w-band, i'-band and 
z'-band are altered to 26.493, 25.641 and 24.843 respec- 
tively. The B-band (F439W) zeropoint is 25.662. 

We have corrected for the small amount of fore- 
ground Galactic extinction toward the HDFN using the 
COBE/DIRBE & IRAS/1SSA dust maps of Schlegel, 
Finkbeiner & Davis (1998). For the GOODS-N field, 

2 available from ftp://archive.stsci.edu/pub/hlsp/goods/ 



Galactic extinction is given by E(B-V) = 0.012 mag, 
equivalent to Apn 5 w — 0.025 and Ap S50 LP — 0.018. 

For consistency, given minor discrepancies between ob- 
ject coordinates in different epochs, the astrometry pre- 
sented in this paper is taken from the world coordinate 
system provided on the third epoch of the GOODS v0.5 
released data. It should be noted that the first epoch of 
data was released with only an approximate astrometry 
solution although this was corrected subsequently. 

Within each epoch of data, object detection was per- 
formed using the SExtractor photometry software pack- 
age (Bertin & Arnouts 1996). Fixed circular apertures 
l'/O in diameter were trained in the z'-band image and 
the identified apertures were used to measure the flux at 
the same spatial location in the i'-band image, running 
SExtractor in 'two-image' mode. This allows for identifi- 
cation of any object securely detected in the z'-band but 
faint or undetected at shorter wavelengths - the expected 
signature of high redshift objects. For object identifica- 
tion, we demanded at least 5 adjacent pixels above a flux 
threshold of 2a per pixel (0.01 counts pixel -1 s _1 ) on the 
drizzled data (with a pixel scale of 0'.'05 pixel -1 ). 

A preliminary catalog of detected objects was cre- 
ated for each epoch of data. In the compilation of the 
catalogues a number of 'figure-eight' shaped optical re- 
flections and diffraction spikes near bright stars were 
masked, as was the gap between the two ACS CCDs and 
its vicinity in each tile. Excluding the masked regions, 
the total survey area in the GOODS-N is 200 arcmin 2 , 
although approximately 20% of this area is only observed 
in half of the available epochs. We estimate that the frac- 
tion of each image contaminated by bright foreground 
objects is small (< 2 per cent). 

To select z ~ 6 galaxies, we adapt the now widely- 
applied Lyman break technique pioneered at z ~ 3 by 
Steidel and co-workers (Steidel et al. 1995, 1996). At 
z ~ 6 the integrated optical depth of the Lyman-a for- 
est is 3> 1 and hence the continuum break in an ob- 
ject's spectrum at the wavelength of Lyman-a is large, 
and can be identified using just two filters: one in the 
Lyman forest region, where there is essentially no flux 
transmitted, and a second longward of the Lyman-a line 
(-Vest = 1216 A). At z ~ 6 Lyman-a emission lies at an 
observed wavelength of A b s ~ 8500 A, in the SDSS-z' 
(F850LP) band, and we use the SDSS-i' (F775W) band 
to cover the Lyman forest wavelengths. As discussed in 
paper 1, the sharp band-pass edges of these SDSS filters 
aid in the clean selection of high-redshift objects (such 
as the z ~ 6 QSOs of Fan et al. 2001), and in paper 2 we 
spectroscopically confirmed an i'-drop selected by this 
technique to be a star-forming galaxy at z = 5.78. 

To quantify this, simulations were performed in which 
an ensemble of galaxies with a distribution of spectral 
slopes and luminosities similar to those of the well- 
studied Lyman-break population at z ~ 3 (Adelberger & 
Steidel 2000) were generated. These were uniformly dis- 
tributed in the redshift range 4 < z < 7 and the effects of 
absorption in the IGM simulated using the Madau (1995) 
prescription. The i' — z' colours of these galaxies were 
then calculated and hence the selection completeness in 
each redshift bin assessed. In these simulations, > 90% 
of galaxies brighter than our magnitude limit at z > 5.6 
and > 99% of galaxies at z > 5.7 satisfied our colour 



HST imaging and Keck Spectroscopy of /-drops 



3 



cut criterion of (i! — z')ab > 1.5. The fraction of z > 4 
interlopers with these colours was negligible(< 1%). The 
luminosity function and distribution in spectral slope of 
galaxies at z ~ 6 are currently unknown (see section T5.4|l 
so these simulations are indicative only but do suggest 
that this technique would successfully recover a popula- 
tion equivalent to that at z ~ 3. 

We impose a magnitude limit corresponding to L* in 
the rest-frame UV for the Lyman-break population at 
z = 3 (Steidel et al. 1999): M* AB {lb00k) = -21.1 mag, 
which translates to z' AB = 25.6 mag at z = 5.8 where 
our selection sensitivity peaks (see section I4.1f) and to 
an unobscured star formation rate of 15 M^yc" 1 (see 
section |SJ. This cut, a la detection in each epoch of the 
GOODS/ ACS data, ensures that each object is securely 
detected in z'. 

Randomly placed artificial galaxies were generated in 
each epoch, using the IRAF. Artdata package, with z' 
magnitudes in the range 20-28 and distributed acording 
to a Schecter luminosity function with a faint end slope 
a = —1.6 and m* = 25.6, the parameters of the Lyman 
break population at z = 3 projected to z — 6. We con- 
firm that at our magnitude cut and surface brightness 
limit (see Figure |3J), we reliably recover >98 per cent of 
such galaxies in the z' band in each epoch, with the re- 
maining 2 per cent lost to crowding (confusion due to 
blended objects). 

In order to minimize the effects of photometric errors, 
un-rejected cosmic rays and transient objects on our final 
candidate list, an independent selection of objects satis- 
fying the criteria of (i 1 — z')xb > 1-5 and z' AB < 25.6 mag 
was made on each epoch of the GOODS data set (see Fig- 
ure^) . This colour cut criterion allows seperation of low 
from high redshift galaxies (subject to photometric scat- 
ter) as shown by figure 2 in paper 1 while minimising 
the contamination due to elliptical galaxies at z ~ 2. In 
order to ensure the completeness of the final sample, ob- 
jects detected with magnitudes and colors within approx- 
imately 2a of our final color and magnitude cut in each 
epoch were also considered. The objects in the resulting 
sub-catalog were then examined for time-varying tran- 
sient behavior across those epochs in which they were 
observed, and a number of objects which satisfied our 
color criterion in individual epochs were rejected at this 
stage. 

The fluxes of the remaining candidate objects were 
then measured directly within a 1" diameter aperture 
in each epoch and these fluxes averaged to give a final 
magnitude, thereby reducing the photometric error for 
each object. To convert to total magnitudes, an aper- 
ture correction of 0.09 mag was applied to all objects, 
determined from aperture photometry of bright but un- 
saturated point sources. Although this correction will 
underestimate the magnitude of very extended sources, 
the candidate objects presented in this paper are all com- 
pact and it has not been neccsary to apply different aper- 
ture corrections to different candidates. 

Although a deeper candidate selection may have been 
possible using coadded images for selection, we have 
found that the advantages of being able to eliminate spu- 
rious and transient sources by epoch-epoch comparison 
are significant. An ideal algorithm for co-adding the im- 
ages would reject all deviant pixels allowing identification 
of fainter objects, however such an algorithm is in prac- 
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Fig. 1. — Color selection of our candidates showing the (i' — 
z')ab > 1.5, z' AB < 25.6 mag selection criteria applied in each 
epo ch and indicating the colors of the candidate objects in table 
12.21 A further two objects lying within la our V — z' cut are 
included for completeness. Note that the catalog presented here 
is not clean and indicates the large number of data artifacts that 
can mimic these i' — z' colours. A number of 'objects' identified 
in each epoch were in fact connected to unmasked satellite tracks 
or optical refraction artifacts and these objects were rejected from 
the final selections. 



tise difficult to obtain with just three or even five epochs 
of data and coadding the data here will lead to the loss of 
time-dependent information or information on poorly re- 
jected spurious pixels. In several cases the time-averaged 
colour of an object was significantly affected by hot pixels 
or poorly rejected cosmic rays in a single epoch. The se- 
lection method outlined above allowed epochs for which 
the data were unreliable to be rejected when making a 
final selection. As a result we have chosen to work at a 
relatively bright magnitude limit on the shallower single 
epoch data, the completeness of our resulting sample is 
discused in section 0] 

In table 12. 2l we present details of the nine objects meet- 
ing our selection criteria. In addition the two objects 
which lie within lcr of being selected are included in the 
following discussions for completeness. Postage stamp 
images of each object in the v, i', z' and HK 1 bands (see 
section l2~3"|) are presented in figureEl None of the i'-drops 
in table !2.2l were formally detected at > 3cr in the F606W 
(v) or F435W {B) band (v AB < 27.3, B AB < 27.0 mag) 
of GOODS/ ACS observations except candidate 6 which 
may be subject to contamination from a nearby object 
(figure 0) . 

As illustrated by figure 3 of the 9 i'-drop objects in 
table l2~2l are unresolved in the GOODS/ ACS z' images 
(R/j =0'.'05). These objects may reasonably be Galactic 
stars, high redshift QSOs or compact galaxies. Previ- 
ous studies have shown that galaxies at very high red- 
shift are barely resolvable with HST (and arc sometimes 
unresolved), examples including the z > 5 candidates 
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of Bremer et al. (2003) which have R h ~0"l-0"3, and 
the barely-resolved spectroscopically confirmed z=5.78 
galaxy SBM03#3 which has an R h =0"08 (Paper 2). 
With this in mind, we consider each object individually 
rather than requiring that all galaxy candidates are re- 
solved. 

2.2. X-ray properties: Non-Detection of i' -Drops by 
Chandra 

To avoid possible contamination of our sample by high 
redshift AGN we examined publically-available deep data 
from the Chandra X-ray telescope. 

All of our candidates lie in the region surveyed in a 2 Ms 
exposure by the Chandra X-ray satellite (Alexander et al. 
2003) and none of these objects were detected by that 
survey, allowing us to place a limit on their X-ray fluxes 
0.5 — 2 keV (soft) and 2 — 8 keV (hard) bands at between 
1.9 x lO" 17 and 9.3 x lO" 17 ergs cm s and between 
1.1 x 10~ 16 to 7.5 x 10~ 16 ergs cm -2 s _1 respectively, 
varying according to the non-uniform exposure across 
the GOODS-N field (D.M. Alexander, private commu- 
nication). A powerful AGN (quasar) is expected to 
have a rest-frame 2-8 keV luminosity > 10 44 ergs -1 
(Barger et al. 2003), which yields a flux at z = 5.8 of 
~ 3 x 10 _16 ergs _1 cm -2 . Thus we would expect to de- 
tect X-ray luminous quasars for most of our sources at 
the relevant search redshifts probed in this paper. 

This renders the identification of any of our candidates 
as powerful high-redshift AGN unlikely although non- 
detection of x-ray flux at these limits does not rule out 
the presence of fainter high redshift AGN such as those 
observed in i=s 3% of Lyman break galaxies at z ~ 3 (Stei- 
del et al. 2002). In this z ~ 3 sample, AGN contribute 
less than 2% of the total rest-frame UV luminosity. Cal- 
culations in section [5] require an implicit assumption is 
that the rest-frame flux is dominated by star formation. 
Although an absence of luminous AGN does not necesar- 
ily imply a UV continuum dominated by star formation, 
it does decrease the likelihood of significantly non-stellar 
origins for the flux. 

2.3. Near-infrared Imaging 

As figure 01 illustrates, a color cut of (i' — z')ab > 1-5 
can be used to select galaxies with z > 5.6. This color-cut 
criterion alone may also select elliptical galaxies at z ~ 2 
(e.g., Cimatti et al. 2002) and cool low-mass M/L/T stars 
(e.g., Kirkpatrick et al. 1999, Hawley et al. 2002). We 
can guard against such contaminants by utilizing near- 
infrared photometry. In particular the relatively blue 
[z! — HK') > 3 color expected of z ~ 2 ellipticals al- 
lows these objects to be distinguished from high redshift 
candidate objects. 

In order to improve our discrimination between high 
and low redshift objects we have supplemented the 
GOODS EST /ACS data with publically-available near- 
infrared images from the Hawaii- HDFN project 3 . The 
HK'-ba,nd data is described in Capak et al. (2003), and 
was obtained on the University of Hawaii 2.2m telescope, 
with a pixel scale of O'.'SpixeP 1 . It covers an area of 
0.11 deg 2 to a depth of HK' AB = 22.5 (encompassing 
the entire GOODS area), with the central 9' x 9' region 

3 available from http://www.ifa.hawaii.edu/~capak/hdf/ 
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Fig. 2.— i'-drops in the GOODS-N field: HST/ACS and Hawaii 
HK' , 4 arcsecond boxes. ACS images are constructed from all 
available epochs. An asterisk indicates that an object forms part 
of our spectroscopic sample. 



imaged to a greater depth of HK' AB — 23.7 (2 a, 1" di- 
ameter aperture). 

SExtractor was again used to identify objects in the 
HK' images with 5 or more adjacent pixels exceeding the 
background noise by 2a. Matches for the i'-drop candi- 
date objects were assigned if positions of the sources were 
within 1", although we carefully examined the z'-band 
image to see if nearby non-z'-drop sources may instead 
be responsible for any near-infrared flux. An aperture 
correction of 0.4 mag was applied to the 1" aperture 
magnitudes in order to estimate the total flux. 

Low redshift z ~ 1 — 2 elliptical galaxies are known 
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ID 


RA 


Declination 


2 AB 




(i' - z')AB 


(z< - HK') AB 




SFR-= 58 


Final 




(J2000) 


(J2000) 










arcsec 


h 7Q 2 MqjV 1 


Cut? 


1* 


12 h 36 m 58?13 


+62°18'51'.'5 


23.05 ± 0.02 


24.96 ± 0.07 


1.91 ± 0.07 


0.21 ± 0.11 


0.051 1 " 


[140] 


No 


2 + 


1 oh i *7m Q /IS f\A 
1Z 61 34.U4 


COOT CQ" A 

+ 62 10 DO. 4 


oq on -L. n no 
Zo.ZU ± U.02 


O/t 70 _L n c\a 
24. (9 ± 0.U6 


i tin _L n ne 
1.60 ± U.06 


1 £IO _L_ n CiA 

1.62 ± 0.U4 


0.053+ 


[122] 


JMO 


3 


12 h 36 m 38?85 


+62°14'51'.'9 


23.91 ±0.04 


25.80 ±0.16 


1.89 ±0.17 


1.89 ±0.05 


0.054t 


[63] 


No 


4 


12 h 35 m 37?21 


+62°12'03'/5 


24.63 ±0.07 


26.13 ±0.25 


1.59 ±0.26 


< 2.22 (2 a) 


0.164 


33 


Yes 


5* 


12 h 36 m 48?47 


+62°19'02'/0 


24.63 ±0.07 


26.29 ±0.27 


1.67 ±0.28 


1.70 ±0.08 


0.197 


33 


Yes 


6* 


12 h 37 m 39?29 


+62°18'40'.'2 


25.06 ±0.10 


26.93 ± 0.54 


1.87 ±0.55 


< 2.65 (2 a) 


0.201 


22 


Yes 


7 


12 h 35 m 52?35 


+62°ll'42'.'l 


25.36 ±0.14 


27.09 ±0.66 


1.73 ±0.67 


< 2.95 (2 a) 


0.083 


17 


Yes 


8 


12 h 36 m 50?41 


+62 o 20'12'.'6 


25.39 ±0.14 


< 27.19 (2a) 


> 1.71 (2 a) 


< 2.98 (2 a) 


0.177 


16 


Yes 


9 


12 h 36 m 48?74 


+62°12'17'.'l 


25.54 ±0.16 


< 27.19 (2a) 


> 1.56 (2 a) 


< 1.93" (2 a) 


0.159 


14 


Yes 


10* 


12 h 36 m 53?79 


+62°11'18'.'2 


24.22 ±0.05 


25.67 ±0.14 


1.45 ±0.15 


1.96 ± 0.07 


0.054+ 


[47] 


No 


11 


12 h 36 m 49?79 


+62°16'24('9 


24.86 ± 0.09 


26.23 ±0.25 


1.37 ±0.26 


< 2.45 (2 a) 


0.058+ 


[26] 


No 



Note. — Objects 10 and 11 (below the horizontal line) lie marginally outside our selection criteria but are included for completeness. 
All magnitudes are measured within a 1 arcsec diameter aperture. The half-light radius, R^, is defined as the radius enclosing half the flux 
in the 2' band. The nominal star formation rate (SFR) is calculated assuming objects are placed in the middle of our effective volume 
(a luminosity-weighted redshift of z = 5.8) and are shown in square brackets if the candidate is not included in our final high redshift 
selection. 

*Keck/DEIMOS spectrum obtained 
t Unresolved point source 

** Additional Near-IR photometry is available for candidate 9 from the deep HST/NICMOS imaging of the HDFN (Thompson et al. 
1999). This object has FIIOWab(J) = 25.10 and F160W A b(#) = 24.60 in this dataset. 
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Fig. 3. — The magnitude— size distribution of our GOODS-N i- 
drops (circles) and objects in the field. Note the cut-off at large 
half-light radii (R/i) due to surface brightness dimming (the solid 
curve is our completeness limit). Our magnitude limit is marked 
with a dotted line. 



to be a contaminant in colour selected samples at z > 5 
(Capak et al. 2003, Bremer et al. 2003). Our i' - z' 
colour cut criterion was designed to exclude these but 
contamination is still possible due to scatter in the pho- 
tometric properties of these objects. Although the near 
infrared colors of many of the objects presented in table 
12.21 as shown on figure 01 are presented only as upper 
limits, those limits all place (V — HK')ab < 2.5 with 
the exception of those for candidates 7 and 8. Given the 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 
(i'-z) colour, AB mags 

Fig. 4. — Evolutionary tracks and color discrimination: the 
solid black line indicates the evolutionary path followed by ellipti- 
cal galaxies with redshift (calculated using the Coleman, Wu and 
Weedman (1980) spectral template). Also indicated by lines are 
evolutionary paths for Sab (dashed), Sbc (dotted) and Im (dot- 
dash) galaxies and for a Kinney et al. starburst model (dot-dot- 
dash). Tick marks on tracks indicate unit redshift intervals be- 
tween 2 = (lower left) and z = 7 (upper right). The vertical dot- 
ted line indicates our basic color- cut criterion of (i' — z')ab > 1-5. 
The colors of candidates in table |2~21 are indicated with error bars 
or as upper limits in (z 1 — HK')\ B . The colors of typical L class 
stars are shown as small diamonds, derived from the spectra from 
Reid et al. (2001). 



red (i' — z')ab colour of candidate 8, only candidate 7 is 
marginally consistent with being a low redshift elliptical 
galaxy. The lower limit of (z' — HK')ab < 2.95 for this 
object, however, does not preclude the possibility that 
candidate 7 may also lie at much higher redshift. 
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The region of colour space occupied by L and M class 
stars is less clearly distinct from that occupied by high 
rcdshift galaxies. Four objects (candidates 2, 3, 5 and 
10) have colors consistent with identification as L class 
stars although all four are also still consistant with a 
high redshift interpretation if the spectral slope of high 
redshift objects is steeper than that observed in z ~ 3 
LBGs. Candidate 1 has colours which are consistant with 
either a high redshift galaxy or a mid-M class star. Near- 
IR colours, together with their spectroscopy and other 
properties, are considered in section 13.21 where a final 
high rcdshift candidate selection is made. 



Table 2. Summary of Spectroscopic Results 



Object Observation 

GOODS-S SBM03#1* Line emission at 8305A - identified as Lyman-a 

GOODS-S SBM03#2 No observed spectral features 

GOODS-S SBM03#3* Line emission at 8245A - identified as Lyman-a 

GOODS-S SBM03#5 No observed spectral features 

GOODS-S SBM03#7* No observed spectral features 

GOODS-N#l Continuum Flux longward of approx. 7200A 

GOODS-N#2 Continuum Flux longward of approx. 7800A 

GOODS-N#5* No observed spectral features 

GOODS-N#6* Line emission at 8804A - identification uncertain 

GOODS-N#10 No observed spectral features 



3. /-DROPS - SPECTROSCOPIC PROPERTIES FOR THE 
HDF NORTH AND CDF SOUTH SAMPLES 

3.1. Spectroscopy - Observation Data 

We have obtained slit-mask spectra of about half our 
i'-drop sample, using the Deep Imaging Multi-Object 
Spectrograph (DEIMOS, Faber et al. 2002; Phillips et al. 
2002) at the Cassegrain focus of the 10-m Keck n Tele- 
scope. DEIMOS has 8 MIT/LL 2fcx4fc CCDs with 15 ^m 
pixels and an angular pixel scale of 0.1185 arcsecpix -1 . 
The seeing was typically in the range 0'.'7 - l'.'O FWHM, 
smaller than or comparable to the slit width of l'.'O. 

The observations were obtained using the Gold 
1200 line mm -1 grating in first order. The grating was 
tilted to place a central wavelength of 8000 A on the de- 
tectors, and produce a dispersion of 0.314 A pixel -1 . We 
sample a wavelength range of approximately AA bs 6600— 
9100 A, corresponding to a rest-frame wavelength of 
AA rest 950 -1300 A at z ~ 6. Wavelength calibration was 
obtained from Ne+Ar+Hg-r-Kr reference arc lamps. The 
spectral resolution was measured to be AAp^HM ~ 1.4 A 
(Awfwhm ~ 55 km s _1 ) from the sky lines. A small 8 A 
region in the middle of the wavelength range for each 
object is unobserved as it falls in the gap between two 
CCDs. The fraction of the spectrum affected by skylines 
is small (< 2 per cent). Flux calibration was determined 
using the spectra of the alignment stars of known broad- 
band photometry (/ sa 16— 19 mag), used to position the 
two slitmasks through 2" x 2" alignment boxes. 

In Paper 2 we already reported the confirmation of a 
z = 5.8 galaxy (SBM03#3 in GOODS-S from Paper 1), 
based on a 5.5-hour DEIMOS spectrum taken on UT 
2003 January 08 & 09. In section l3~3l we present spectra 
from other GOODS-S i'-drops on the same slitmask. Our 
spectroscopy reveals a second z = 5.8 Lyman-a emitter 
in the same field. The observations and data reduction 
for this slitmask are as detailed in paper 2. 

In addition, we have recently obtained ultra-deep spec- 
troscopy of our GOODS-N i'-drops over the 5 nights of 
U.T. 2003 April 2-6. We used several slitmasks to target 
a subset of five of the eleven objects in table 12.21 (can- 
didates 1, 2, 5, 6 and 10). The i'-drops were placed on 
slitmasks for a primary program targetting z ~ 1 E/S0 
galaxies in the GOODS-N field (Treu et al. in prep). The 
selection of i'-drops from table 12.21 for the masks was 
purely geometric and accordingly randomized. The po- 
sition angle of slits on the mask was chosen to be 45.3° 
so as to align the large (16.5 x 5arcmin) DEIMOS field 
of view with the axis of the GOODS-N EST /ACS obser- 
vations. The spectra were pipeline reduced using version 



Note. — Objects identified as high redshift candidates are 
marked with an asterisk. 

1.1.3 of the DEEP2 DEIMOS data reduction software 4 . 

A total of 37.8 ksec (10.5-hours) of on-source integra- 
tion was obtained in each of three slitmasks, and this 
was broken into individual exposures each of duration 
1800 s. These spectra are therefore a factor of V2 deeper 
than the 5.5-hour spectra of GOODS-S introduced in pa- 
per 2. In 10.5 hours of DEIMOS spectroscopy we reach a 
flux limit of 2 x 10 -18 ergs cm -1 s -1 (5 a) for an emission 
line uncontaminated by sky lines and extracted over 8A 
(300 km s -1 ) and 1". 

The results of this spectroscopy are summarised in ta- 
ble [5] and discussed in more detail in the following sec- 
tions. 

3.2. GOODS-North Candidate Spectra and Final 
Selection 

Since all the objects presented in table 12.21 satisfy our 
primary i' — z' color selection criterion and do not have 
z' — HK' colors indicative of being lower redshift ellipti- 
cal galaxies, our most significant remaining source of con- 
tamination in a high redshift sample is likely to be cool 
galactic stars. These constitute two of the nine z-drop 
objects in the GOODS-S field (paper 1) and are often 
excluded from colour selected samples by the simple ex- 
pedient of rejecting all unresolved sources (e.g. Bouwens 
et al. 2003). In this section we consider the interpretation 
of the candidate spectra for those objects for which such 
is available, the probable classification of each candidate 
in table 12.21 and the construction of a robust sample of 
high rcdshift galaxies. 

Five objects in GOODS-N were targeted for spectro- 
scopic observations. Two yielded continuum detections, 
and a likely identification as cold stars (candidates 1 and 
2). One spectrum (candidate 6) shows a single extended 
emission line feature at 8804 A, with the peak emission 
spatially offset from the continuum location by ~ 2" (and 
hence this spectral feature may be associated with an ad- 
jacent foreground object). No signal was detected in the 
other two spectra (candidates 5 and 10). A description 
of each spectrum follows. 

3.2.1. GOODS-N i> -Drop #1 

The spectrum of GOODS-N i'-drop #1 does not ap- 
pear to show any distinctive emission lines. The object 

4 available from http : //astron.berkeley . edu/~cooper/deep/spec2d/ 
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continuum, however, is clearly detected at long wave- 
lengths. Figure presents a box-car smoothed spectrum 
of this object together with the spectrum of an L-class 
star (2MJ0345432+254023, class LO for comparison. 

The object's spectrum and compactness =0'/05 (un- 
resolved) are both strongly suggestive that this object 
may be identified as an M or early L class star and al- 
though the spectrum obtained has very low signal-to- 
noise, there are suggestive commonalities with the LO 
spectrum such as the absorption feature at ~8400A. In 
addition its colors are marginally consistent with identi- 
fication as a mid M-class star, although the near-infrared 
colors are affected by a near neighbor. We provisionally 
identify this unresolved source as a galactic star, and 
hence exclude it from our final high redshift sample. 

3.2.2. GOODS-N i'-Drop #2 

As in the case of GOODS-N i'-drop #1, the spectrum 
of GOODS-N i'-drop #2 does not appear to show any 
distinctive emission lines but the object continuum is 
clearly detected at long wavelengths. Figure presents a 
smoothed spectrum of this object. Both the continuum 
colors of this object and its compactness (Rh =f/.'05, un- 
resolved) suggest that this object may be an L dwarf 
star, although it lacks the absorption at ~8400A often 
seen in L class stars. The low signal-to-noise spectrum 
of this object shows a moderate cross correlation with 
those of L class stars and a weaker correlation with the 
composite LBG spectrum of Shapley et al. (2003), red- 
shifted to z ~ 6. We tentatively identify this object as 
a cool galactic star and it is therefore excluded from our 
final list of good high redshift candidates. As the classifi- 
cation of this object is uncertain it should be noted that, 
as one of the brighter objects selected by our (£' — z')ab 
color criterion, this object would have a significant effect 
on our calculations in section [5] if it was shown to lie at 
high redshift, as is discussed in that section. 

3.2.3. GOODS-N i'-Drop #5 

Although a spectrum was obtained for GOODS-N i'- 
drop #5, no spectral features or continuum flux were 
observed. However this object is relatively faint (z' AB = 
24.6) and so we would be unlikely to detect continuum 
emission. Our limit on the equivalent width of an emis- 
sion line is W bs < 10 A (5 u, 300 km s -1 width), provided 
it lies away from a sky line. GOODS-N i'-drop #5 is rel- 
atively extended and so unlikely to be a galactic star. Its 
near-infrared color is also inconsistent with identification 
as an elliptical galaxy at z ~ 2 although the colors are 
not fully consistent with a high redshift identification ei- 
ther. We include candidate 5 in our list of high redshift 
candidates although we note the ambiguity in its near- 
infrared colors. 

3.2.4. GOODS-N i'-Drop #6 

The 2D spectrum of GOODS-N i'-drop #6 showed one 
significant emission line detected at « 20 a at a central 
wavelength of A obs ~ 8804 ± 1 A. Both 2D and ID ex- 
tracted spectra for this line are shown in figure El cen- 
tered on the spatial location of the i'-drop object and the 
central wavelength of the emission line. The emission line 
is spatially extended, stretching from the location of the 
i'-drop in the ACS imaging to peak in intensity approx- 
imately 2" away . 




7000 7500 8000 8500 
observed wavelength / A 



9000 



Fig. 5.— GOODS-N i'-drops #1 and #2 — Spectra are binned 
into 10 A bins (33 pixels) and are extracted over 1" (9 pixels). The 
spectrum of candidate 2 is offset by 2 X ergs cm" 2 s _1 for 

clarity. The spectrum of the L0 Dwarf star 2MJ0345432+254023 
(2MASS, Reid et al. 1999) is plotted for comparison (offest by 
1 X 10 — 18 ergs cm -2 s — 1 ), along with the sky spectrum (both with 
arbitrary scaling). The primary spectral features of L class stars 
are marked with arrows and correspond to the following lines re- 
spectively: TiO, TiO, Ki, Ki, Rbi, Rbi, Nai, Nai, TiO, TiO, CS I 
(Kirkpatrick et al. 1999). 



The line flux is (7 ± 2) x 10 18 ergs cm 2 s 1 measured 
between the zero power points at 8798. OA and 8813. 5A 
in an extraction width of 3" centered on the peak of the 
line emission. As figureElshows, the line profile is clearly 
resolved and shows velocity structure, although it is not 
obviously a doublet. The FWHM of the emission line 
is (10.0 ± 1.0)A although the central peak of emission is 
narrower with a FWHM of (4.0 ± 0.25)A. 

GOODS-N i'-drop #6 is a relatively extended, faint 
galaxy with a brighter neighbor at approximately 2" sep- 
aration along the slit direction. The neighboring object 
is expected to lie just off the slit, although given the see- 
ing conditions of this observation (~ 0'.'8) we may expect 
some light from this object to fall into our slit near the 
location of the emission line. As a result we suggest two 
possible interpretations of this observation. 

One possible scenario is that the line emission is asso- 
ciated with this lower-redshift (non-i'-drop) neighboring 
galaxy. The observed line is isolated, with no other lines 
detected. The neighbor has a photometric redshift of 
z = 2.605 ± 0.2, based on its colors in the GOODS/ ACS 
imaging and determined using the template fitting rou- 
tine l HyperZ' (Bolzonella et al. 2000). The observed 
line cannot be easily identified with any emission line 
at z = 2.6 (A rcs t ~ 2442A). It is inconsistent with iden- 
tification as a z = 1.36 [On] (A rcst = 3726,3729A), and 
no doublet structure is obvious, although the line could 
conceivably be Mgn (A rost ~ 2795A) at z = 2.15. 

The second scenario is that we are seeing spatially- 
extended Lyman-a emission at z = 6.24, as has been 
observed in a number of star forming galaxies. Ob- 
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Fig. 6.— (a) GOODS-N i'-drop #6 2D spectrum — Candidate 
Lyman-Q emission line at z = 6.24 (lower panel), sky spectrum 
(upper panel), 2D spectrum boxcar smoothed over 3 pixel resolu- 
tion elements (centre panel). The 2D spectra span a spatial range 
of ±2.5" around the i'-drop location on slit. The emission line is 
centered l'/5 from the candidate object but extends to its location. 
(Right) HST/ACS image showing the slit orientation. Box length 
is 10". As can be seen the neighbouring object lies just off the slit 
but some flux from this neighbour would be expected to fall down 
the slit in ground-based seeing, (b) The ID spectrum extracted 
over 2" around the emission line center, and binned into indepen- 
dent spectral resolution elements of 1.2 A (4 pixels). The features 
around 8790 A and 8830 A are sky-subtraction residuals. 



servations of a lensed z — 4 galaxy by Bunker et al. 
(2000) identified Ly-a emission offset by f-2" from the 
rest frame UV continuum region (although this corre- 
sponds to only 0.7 h^ kpc at z = 4 given the lensing 
amplification of ~ 10 for this galaxy). Similarly Stei- 
del and coworkers (2000) identify two Lyman-a 'blobs' 
associated with but not centered on continuum emission 
Lyman break galaxies at z ~ 3. These Lyman-a regions 
extend over > 100 kpc and are similar to extended 
Lyman-a emission seen in the vicinity of high redshift 
radio-galaxies. 

At z ~ 6 a spatial offset of 2 " corresponds to a physi- 
cal separation of 11.4 h^Q kpc, larger than that observed 
in the z — 4 galaxy but much smaller than the extended 
Lyman-a observed by Steidel et al. (2000). If we identify 
the emission line as Lyman-a then the emission wave- 
length is A res t ~ 1216A giving a redshift for candidate 
6 of z = 6.24. We tentatively adopt this identification 
and candidate 6 remains part of our high redshift galaxy 
sample. 

3.2.5. GOODS-N %' -Drop #10 




Fig. 7.— GOODS-N i'-drop #9 in the ACS z'-band image and 
in the NICMOS FHOWab(J) and F160W A b(#) imaging of the 
Hubble Deep Field North. 



GOODS-N z'-drop #10 was selected with marginal 
i' — z' colors, slightly below our strict color cut criterion. 
The spectra of this object shows no significant features 
or continuum flux in its spectrum. Both its compact 
(unresolved) Rh and (z' — HK')ab color suggests that 
this may be a low mass star. Given that this object did 
not meet our strict criteria, the relative brightness and 
its lack of strong high redshift emission features we ex- 
clude this object from our final selection as a probable 
low mass star although, again, we consider the effects of 
including this object in section [S] 

3.2.6. GOODS-N i' -Drops for Which No Spectra Are 

Available 

In common with GOODS-N i'-drops #2 and #10, 
GOODS-N i'-drop #3 is unresolved and its z' - HK' 
color is completely consistent with identification as an 
L-class dwarf star (Reid et al. 2001). It is thus excluded 
from the final robust high redshift sample although the 
effects of including all our stellar candidates as possible 
high redshift objects is discussed in section [SJ 

GOODS-N i'-drops #4 and #8, on the other hand, are 
relatively extended objects and so are inconsistent with 
identification as galactic stars. Their near-infrared colors 
are also inconsistent with elliptical galaxies at z ~ 2 sug- 
gesting that they may safely be identified as high redshift 
candidates. These candidates remain in our selection. 

GOODS-N z'-drop #9 is relatively extended and can 
be identified with Fernandez-Soto, Lanzetta and Yahil 
(1999) object 311 which was given a photometric red- 
shift of Zphot = 5.64 by those authors, consistent with 
our result. This is the only high-rcdshift candidate object 
here to appear in the Fernandez-Soto, Lanzetta and Yahil 
(1999) photometric redshift catalogue for the HDFN. 
Since GOODS-N i'-drop #9 is inside the central Hub- 
ble Deep Field North region, it is observed in deep NIC- 
MOS near-infrared imaging (Thompson et al. 1999), as 
shown in figured from which we have obtained this ob- 
ject's near infrared magnitudes in the F110Wab(</) and 
F160W A b(#) bands. Its colors of [z 1 - J)ab ~ 0.53 
and (J — H)ab ~ 0.50 are consistent with a high red- 
shift interpretation and inconsistent with the colors of 
low redshift galaxies. As a result we include it in our list 
of final candidates. 

Since GOODS-N i'-drop #11 was selected with 
marginal (i' — z')ab colors, not meeting the strict (i' — 
z')ab < 1-5 color cut, and is also unresolved it has also 
been rejected as a probable low mass star despite the 
lack of a strong near- infrared color constraint. This is 
the only object excluded from our final candidate selec- 
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tion primarily on grounds of its unresolved nature. As 
such we note that its inclusion in our final selection would 
have increased the observed star formation rate of that 
sample by 18% (see section |3J). 

3.3. Spectroscopy of i' -drops in GOODS-South Field 
3.3.1. GOODS-S SBM03#1 

Our 5.5-hour spectrum of GOODS-S i'-drop #1 in pa- 
per 1 (hereafter SBM03#1) reveals a single emission line 
at 8305.3 ± 0.9 A at the same location on the l"-wide 
longslit as the i'-drop (a 20 oo = 03 ,1 32 m 40.0 s , 5 2 ooo = 
-27°48'15.0"). The emission line falls betwe en two sky 
lines, and the 2D spectrum is shown in figure l3~3~T1 

The spectrum is very similar to that of GOODS-S 
SBM03#3 reported in Paper 2, where we argued that 
the most likely interpretation of the solo emission line at 
8245 A was Lyman-a at z = 5.78, given the t'-drop se- 
lection and the fact that we would resolve the [Oil] dou- 
blet at lower redshift (and most other prominent emis- 
sion lines would have neighboring lines within our wave- 
length range). Applying the same logic to our spectrum 
of GOODS-S SBM03#1, it seems likely that the emis- 
sion line is Lyman-a at z = 5.83 - a view supported by 
the characteristic asymmetry in the emission line spec- 
trum, with a sharp blue- wing cut-off 5 . The redshift of 
SBM03#1 is very close to the redshift of SBM03#3 (the 
line centers are separated by 2000 kms -1 ), and the spa- 
tial separation of 8 arcmin corresponds to a projected 
separation of only 2.9 hj^ Mpc. 

The flux in this emission line is comparable to that 
of SBM03#3 at z = 5.78, and is within 30 per cent 
of / = 2 x 10~ 17 erg cm -2 s _1 , extracting over 17pixels 
(2 arcsec) and measuring between the zero-power points 
(8300 - 8318 A). This is a lower limit due to potential 
slit losses, although these are likely to be small, since 
like SBM03#3, SBM03#1 is compact (R u = 0.1", pa- 
per 1). After deconvolving with the ACS z'-band point 
spread function, the half-light radius corresponds to a 
physical size of only 0.5 /i^q 1 kpc. 

In continuum, SBM03#1 is half the brightness of 
SBM03# (paper 2) from the z'-band magnitude. The 

equivalent width is W^J S " — 30±10 A, at the upper end of 
the distribution observed in z = 3—4 Lyman break galax- 
ies (Steidel et al. 1999). The star formation rate from the 
rest-frame UV continuum is 17 Mq yr _1 (paper 1). 

3.3.2. Further Candidates in GOODS-S 

As described in Paper 2, DEIMOS spectroscopy of 
three further candidates (SBM03#2,5,7) was obtained 
with the Keck II telescope. Inspection of these spectra 
reveals no obvious Lyman a emission or continuum emis- 
sion within the wavelength range sampled. At present 
the nature of these sources is unfortunately unclear. 

4. HIGH REDSHIFT GALAXIES - DERIVED PARAMETERS 
AND INTERPRETATION 

4.1. Survey Volume and Completeness 

Galaxies in the redshift range 5.6 < z < 7.0 are se- 
lected by our color-cut provided they are sufficiently lu- 

5 This object has been independently confirmed as lying at 
z = 5.8 by Dickinson et al. (2003) and corresponds to their ob- 
ject SiD002. 



minous. However, as we discuss in paper 1, our sen- 
sitivity to star formation is not uniform over our sur- 
vey volume: at higher redshifts, we are only sensitive to 
more luminous galaxies at our apparent magnitude limit 
of z' AB — 25.6, exacerbated by the increasing fraction of 
the z'-band extincted by the Lyman-a forest, and there 
is a fc-correction as we sample progressively shorter rest- 
frame wavelengths. As a result of this luminosity bias, 
the average redshift of objects satisfying our selection 
criteria is towards the lower end of the redshift range 
above and we calculate this luminosity-weighted central 
redshift as z — 5.8 under the assumption of no evolution 
in the luminosity function from that at z ~ 3. In paper 
1 we followed the methodology of Steidel et al. (1999) to 
quantify the incompleteness due to luminosity bias with 
redshift by calculating the effective volume of the survey. 
If there is no change in L* LBG (1500 A) or a from z = 3 
to z 6 (i.e. no evolution in the galaxy luminosity func- 
tion from that of Steidel et al. 1999 at z — 3), we showed 
that for our ACS z'-band limit for the GOODS data, the 
effective volume is 40 per cent that of the total volume 
in the interval 5.6 < z < 7.0. For our 200 arcmin 2 survey 
area (including areas only surveyed in half of the epochs 
and after excluding the unreliable regions close to the 
CCD chip gaps) this corresponds to an effective comov- 
ing volume of V c s = 2.5 x 10 5 h^Q Mpc 3 . This assumes a 
spectrum flat in /„ at wavelengths longward of Lyman-a 
(i.e., fx oc A _/3 where (3 — 2.0). For a redder spectrum 
with (3 — 1.1 (the mean reddening of the z ss 3 Lyman 
break galaxies, Meurer et al. 1997) the effective volume 
is 36 per cent of the total (V eS = 2.3 x 10 5 /if 3 Mpc 3 co- 
moving) - which would result in an 8% rise in the mea- 
sured star formation density (section [5j compared with 
/3 = 2.0. 

The above discussion, however does not take into ac- 
count the effects of survey incompleteness due to photo- 
metric scatter, extreme intrinsic color variation or cos- 
mological surface brightness dimming. 

Our selection procedure is robust to the effects of pho- 
tometric scatter. By utilizing four epochs of data and 
selecting objects in each which satisfy our criteria, only 
those objects lying more than la from their true mag- 
nitudes in four independent measurements will be en- 
tirely missed (fewer than 1% of the population assum- 
ing Gaussian errors on the magnitudes). It is possible, 
however, that the effects of photometric errors in one or 
more epochs may sufficiently affect the combined flux of 
the observations to take an object below our (i 1 — z')ab 
selection limit. In practise only two objects in our epoch- 
by-epoch selection sub-catalogue showed (i'—z') ab colors 
that consistantly exceeded (i' — z')ab > 1-3 and fluctu- 
ated around our color cut and these were included in 
table as candidates 10 and 11. Both these objects 
were excluded from the final high redshift candidate se- 
lection on the basis of their stellar colors and full-width 
half-maxima. Since relaxing our colour cut criterion to 
(i' — z')ab > 1-3 does not significantly add to the number 
of candidate objects, while at the same time increasing 
the likely fraction of low redshift contaminants (see fig- 
ure we estimate that the fraction of objects falling 
out of our selection due to photometric errors is small. 

The effect of extreme scatter in the distribution of in- 
trinsic object colors on the completeness of our selection 
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Fig. 8. — (a) The lower panel shows the two-dimensional spectrum of GOODS-S SBM03#1 (candidate 1 in paper 1) around Lyman-a 
(100 A by 5arcsec), the middle panel shows this smoothed with a Gaussian of a = 1 pixel, and the upper panel shows the sky spectrum 
for this range. The line is centered at 8305 A. (b) The one-dimensional spectrum of SBM03#1 around Lyman-a, extracted over a 9pixel 
(laecsec) width. The data have been binned into independent spectral resolution elements of 1.3 A (4 pixels). 



is difficult to quantify given the unknown nature of the 
z ~ 6 population. Our color selection criterion illus- 
trated in figure 0] is based upon the colors predicted by 
the evolution of empirical galaxy templates. This pro- 
cedure has been shown to work well in the selection of 
Lyman break galaxies at z ~ 3 — 4 but does not account 
for atypical galaxies or for any significant deviation in 
the spectral shape of very high redshift objects from the 
templates. 

As described in section simulations were made to 
assess the completeness of our i' — z' colour selection cri- 
terion on a galaxy population with the characteristics 
of the z ~ 3 Lyman break galaxies (Adelberger & Stei- 
del 2000). These found that > 99% of z > 5.7 galaxies 
brighter than our luminosity limit in such a population 
would be identified using this technique, with a some- 
what lower completeness in the 5.6 < z < 5.7 redshift 
bin where the completeness depends sensitively on the 
filter transmission and the model of IGM absorption. 
Relaxing the colour cut criterion to (i' — z')ab > 1-3 
would increase our completeness in the 5.6 < z < 5.7 
redshift bin by ss 10% without affecting completeness in 
the higher redshift bins but, as noted above, would also 
slightly increase the levels of contamination from z ~ 2 
elliptical galaxies. 

The characteristics of the z ~ 6 population are cur- 
rently unknown with neither the luminosity function or 
the distribution of spectral slopes accessible to the avail- 
able data and it is possible that there is significant evo- 
lution in these properties in the interval 3 < z < 6. How- 
ever, the %' — z 1 colour of galaxies at these high redshifts 
is driven primarily by the break caused by absorption 
in the intervening intergalactic medium. As a result the 
i 1 — z' colour is effectively independent of spectral slope 
and the colours of a wide variety of intrinsic spectral en- 
ergy distributions converge for z > 5. 

The single object confirmed to lie within our redshift 
range in the GOODS-N field (HDF 4-473 at z = 5.60, 
Weymann et al. 1998) is slightly too faint to appear in 
our object selection {z' AB = 26.66, Bouwens et al. 2003a). 
More significantly, this galaxy is also too blue to ap- 
pear in our selection ((i' — z')ab = 1-2, Bouwens et al. 
2003a), having been observed as a V band drop-out, and 



so raises the possibility that we are excluding real objects 
by applying our color-cut. We note that HDF 4-473 is 
at the very lower end of the redshift range to which we 
are sensitive, and hence in a bin where our completeness 
is relatively poor. Nonetheless this object highlights the 
caution which must be applied to the interpretation of 
any sample selected purely on photometric criteria. 

4.2. Physical Sizes of the i' -drops and Surface 
Brightness Selection Effects 

Our final catalog selection includes objects with ob- 
served sizes in the range R/, ~ 0'.'08 — 0'.'20, similar to 
the observed range of sizes in the 6 z ~ 6 candidates 
we presented in paper 1 . This range corresponds to pro- 
jected physical sizes of 0.3 — 1.0/i^j kpc at z ~ 5.8, de- 
convolved with the instrument PSF. The size range of 
our total sample of 12 z ~ 6 objects (paper 1 + this 
work) is comparable to, or smaller than, those of z > 5 
galaxies reported by Bremer et al. (2003). Bremer and 
coworkers identified a sample of 44 sources with half-light 
radii of R/j ~ O'.'l — 0'.'3 and suggested that, when taken 
in conjunction with observations at lower redshifts, their 
results implied a modest decrease in galaxy scale lengths 
with redshift. Our results appear to support this conclu- 
sion, although it should be noted that at = 25.6, 
objects more extended than R^ ~ 0'.'25 (1.5 h^Q kpc) 
fall out of our selection due to low surface brightness 
(as shown in figure (5J, thus placing a significant obser- 
vational selection effect on the observed distribution of 
sizes. 

Simulations were performed using the IRAF. Artdata 
package, placing artificial galaxies in the GOODS im- 
ages in order to test the recoverability of such objects 
given our SExtractor detection parameters. We model 
the input population as having exponential disk surface 
brightness profiles with half light radii distributed uni- 
formly between the instrumental point spread function 
(0705) and 075. We find that we reliably recover >95% 
of the input galaxies for half-light radii Kh < 0725 and 
>98% of galaxies with R^ ~ 0720 with the remaining 2% 
being lost because of crowding. 

Nonetheless, at brighter magnitudes our objects do not 
occupy the full range of half-light radii allowed by the ob- 
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servational constraint: at z' AB — 24.7 (the magnitude of 
our brightest high redshift candidate) we are sensitive to 
objects as large as R/j ~ 0'.'34 (see figure This sug- 
gests that most objects at high redshift are small enough 
to fall within our surface brightness selection limit. 

Star forming regions in the local universe are often very 
compact in the UV continuum and it is likely that this 
is also true at high redshift. Lowenthal et al. (1997) re- 
port that Lyman break galaxies at z ~ 3 are compact 
and, if projected to z ~ 3, this population would range 
in half-light radius from 0'.'08-0'.'35 with a median size of 
Rh = 0'.'18. Although just over half of this range of half- 
light radii would be detectable at our magnitude limit, 
the full range of half-light radii are accessible to this sam- 
ple at brighter magnitudes and failure to detect such ob- 
jects may indicate a reduction in the physical scale of 
galaxies with increasing redshift. At z > 5 Bremer and 
coworkers (2003) identified only 2 candidate objects on 
scales significantly larger than R/j ~ 0'.'3 (1.8 kpc at 
z ~ 5.2) in a sample of 44 galaxies and considered these 
objects most likely to be z ~ 1 red galaxies. Iwata et al. 
(2003) made a similar selection of Lyman break galax- 
ies at z > 5 and found 40 per cent of their candidate 
objects were unresolved although they were limited by 
seeing conditions in their ground based observations. As 
Bremer et al. discuss, they expect a contamination from 
lower redshift galaxies of ~ 20% in their sample and also 
expect such objects to be resolved on scales of > 0'.'3. 

While this constraint on the physical size of observed 
objects is not explicitly accounted for in our effective vol- 
ume calculation, the change in the angular diameter of 
objects of the same physical size over our redshift range 
is only 10 per cent. Given the evidence for a slight reduc- 
tion in intrinsic object size with redshift (Bremer et al. 
2003, Roche et al. 1998), and the failure of our candidates 
to occupy the full size range allowed by our observations 
even at fainter magnitudes, we expect the fraction of ob- 
jects lost due to surface brightness effects to be small 
(< 2%) and the effect on our effective survey volume to 
be negligible. 

5. THE STAR FORMATION RATE AT Z ~ 6 
5.1. The Global Star Formation History in GOODS-N 

Since the primary sources of rest-frame ultraviolet pho- 
tons within a galaxy are short-lived massive stars, UV 
flux is a tracer of star formation and its rest-frame UV 
flux density can be used to quantify the total star for- 
mation rate (SFR) within a galaxy. This rest-frame ul- 
traviolet continuum is measured at z ~ 6 by the z' pho- 
tometric band. 

The relation between the flux density in the rest-UV 
around ~ 1500 A and the star formation rate (SFR in 
M yr~ 1 ) is usually assumed to be given by Lyv = 
8 x 10 27 SFRergs s" 1 Hz" 1 (Madau, Pozzetti & Dickinson 
1998) for a Salpeter (1955) stellar initial mass function 
(IMF) with 0.1 M Q < M* < 125 Mq. 

Our limiting star formation rate as a function of red- 
shift was considered in detail in paper 1, section 4 and 
illustrated in figure 6 in that paper. Accounting for filter 
transmission, the effects of the intervening Lyman-a for- 
est and for small fc-corrections to A res t = 1500 A from 
the observed rest-wavelengths our limiting magnitude 
{z' AB = 25.6) we should detect unobscured star formation 



rates as low as 15 [16.5] ft. 70 2 M©yr 1 at 5.6 < z < 5.8 
and 21 [25] Mq yr" 1 at z = 6.1 for spectral slope 
(3 = —2.0 [—1.1], appropriate for an unobscured starburst 
and a redder slope appropriate for mean reddening of the 
z«3 Lyman break galaxies (Meurer et al. 1997) respec- 
tively. 

Steidel et al. (1999) used a spectroscopic sample of Ly- 
man break galaxies (LBGs) to constrain the luminosity 
function of these objects around A rcs t = 1500 A obtaining 
a characteristic magnitude at the knee of the luminosity 
function of m* R — 24.48 at (z) = 3.04, with a faint end 
slope a — —1.6 and normalization cj>* « 0.005 h^ Mpc -3 . 
A second sample of galaxies at (z) = 4.13 was consis- 
tant with no evolution in this luminosity function. If 
we further assume that the luminosity function does not 
evolve between z ~ 3 and z ~ 6 then our catalogue 
limit at z' AB < 25.6 would include galaxies down to L^ BG 
(corresponding to a characteristic star formation rate of 
SFR T *j V = 15/ifo Moyr" 1 ) at z ~ 6 (although see sec- 
tion for a discussion of this assumption). 

In the ninth column of table 12.21 we use the ultravi- 
olet flux — star formation rate relation assumed above 
to estimate the inferred star formation rate for each of 
our i'-drop candidates assuming our luminosity weighted 
central redshift (z = 5.8) for each object. In sec- 
tion however, we commented on the nature of each 
of the i'-drop objects and selected a subsample of six 
candidate high-redshift star-forming galaxies with z > 
5.6. These have an integrated star formation rate of 
(135 ± 55) M Q yr" 1 (where the error is based purely 
on Poisson statistics). 

If the luminosity function of Lyman break galaxies is 
assumed not to evolve between z = 3 and z ~ 6 this 
gives a comoving volume-averaged star formation den- 
sity between z = 5.6 and z = 6.1 of psfr = (5.4 ± 2.2) x 
10~ 4 h 70 Mq yr -1 Mpc -3 for objects with star formation 
rates > 15 Mq yr -1 . The inclusion of all objects in 
table 1 (the six objects in our robust high redshift sam- 
ple and the five probable stars) would give a value some 
4x higher than this, with half of the integrated star for- 
mation rate contributed by the unresolved low mass star 
candidates #1 and #2 (see section l3~2*|) . While it is true 
that none of our five unresolved candidates can be con- 
clusively identified as galactic stars, we are confident that 
this figure would seriously over-estimate the true value 
given the demonstrated stellar contamination of similar 
samples (e.g. paper 1). 

Star formation densities calculated for both our high 
redshift gala xy selection (filled square) and for all i-drops 
in table 12.21 (open square) are compared to global star- 
formation rates found in paper 1 and at lower redshifts 
in figure El adapted from Steidel et al. (1999) and re- 
calculated for our A cosmology and higher limiting star 
formation rate. If our high redshift candidate selection is 
adopted, then for Lyman break galaxies with star forma- 
tion rates > 15 Mq yr -1 (L*_ 3 ) it appears that the 
observed comoving star formation density was ~ 9 times 
less at z ~ 6 than at z = 3 (based on the bright end of the 
UV luminosity function), a conclusion reinforced by the 
close agreement between our results for the GOODS-N 
and GOODS-S data. If all i-drops in table !2.2l are consid- 
ered the fall in star formation density is less marked (a 
factor of 2), however given the significant contribution 
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the five unresolved objects make to the star formation 
rate and the evidence against their identification as high 
rcdshift objects, we are reluctant to bias our results un- 
duly high by their inclusion in the star formation density 
estimate presented here. 

The close proximity, both in velocity and in projected 
space, of the two confirmed z ~ 5.8 galaxies is suggestive 
of the presence of a group of galaxies at this redshift in 
the GOODS-S field. Nonetheless, the similar results ob- 
tained in both fields suggests that our survey fields are 
sufficiently large that cosmic variance, a problem that 
has always plagued small deep- field surveys, does not 
contribute significantly towards the uncertainties associ- 
ated with the determination. 

The observation that there is evolution in the lumi- 
nosity function of Lyman break galaxies in the redshift 
interval 3 < z < 6 (as observed by the bright end, 
SFR> 15 /iy 2 Mq yr _1 ) is supported by simulation. If 
an ensemble of galaxies with properties identical to the 
z ~ 3 Lyman break population in distribution of spec- 
tral slope and luminosity function (Adelberger & Steidel 
2000) is projected to lie at z > 4 then we would expect to 
detect 49 ± 10 such galaxies at z > 5.6 to our magnitude 
limit, where the error bars incorporate contributions of 
similar magnitude from poisson sampling noise and from 
uncertainties in the z ~ 3 luminosity function. Hence 
the observed fall in number density of bright galaxies is 
strongly suggestive of evolution in either the normalisa- 
tion ((/>* , and hence star formation rate) or in the shape 
(L* or a) of the Lyman break luminosity function in this 
interval (see section IB~4l for further discussion). 

Previous authors have integrated the star forma- 
tion density over the luminosity function for their ob- 
jects down to star formation rates of 0.1 L* (SFR= 
1.5 h^Q Mq yr _1 ) even if these SFRs have not been ob- 
served. We prefer to plot what we actually observe (SFR 
> 15 Mq yr" 1 ) and not to extrapolate below the 
limit of our observations. Unfortunately, calculation of 
the effective survey volume requires an implicit assump- 
tion of the shape of the very uncertain z ~ 6 Lyman 
break galaxy luminosity function, as discussed in section 
14.11 If the Steidel et al. (1999) luminosity function for 
z ~ 3 LBGs (used to calculate this volume) is assumed to 
hold at z ~ 6 and is integrated down to 0.1 L* z=3 then the 
global SFR is increased by a factor of ~ 5 for a = —1.6. 
See section I5~4l for a discussion of the effects of evolution 
in the luminosity function of these objects. 

5.2. Extinction Corrections 

The rest-frame UV is, of course, very susceptible to 
extinction by dust, and such large and uncertain correc- 
tions have hampered measurements of the evolution of 
the global star formation rate (e.g., Madau et al. 1996). 
Although extinction corrections increase the UV flux, 
and hence the implied star formation rate, the method 
used here to select galaxies (the Lyman break technique) 
has been consistently applied at all redshifts above z w 2. 
Thus equivalent populations should be observed at all 
redshifts and the resulting dust correction alters the nor- 
malization of all the data points in the same way. Hence 
the fall in the volume-averaged star formation rate be- 
tween z ~ 6 and z rs 3 derived from our i'-drops does 
not depend on dust extinction, unless the dust content 
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Fig. 9. — The 'Madau-Lilly' plot illustrating the evolution of 
the comoving volume-averaged star formation rate. Our work is 
plotted as solid symbols: the result from paper 1 (GOODS-S) 
is shown as a filled circle and from this work (GOODS-N) as a 
filled square (note that this point includes only objects in our final 
high redshift selection). The large open square at z = 5.8 indi- 
cates the integrated star formation rate in the GOODS-N if the 
five unresolved galactic star candidates are assumed to be unre- 
solved galaxies at z ~ 6. Other determinations have been recal- 
culated for our cosmology and higher star formation rate limit of 
15h~Q Mq yr" 1 (instead of 1.5M yr- 1 for Steidel et al. 1999), 
assuming a slope of a = —1.6 for z > 2 and a = —1.3 for z < 2. 
Data from the CFRS survey of Lilly et al. (1996) are shown as 
open circles; data from Connolly et al. (1997) are squares; and 
the Lyman break galaxy work of Steidel et al. (1999) is plotted as 
crosses, of Bouwens et al. (2003b) by pentagons, of Fontana et al. 
(2002) as inverted triangles and that by Iwata et al. (2003) as an 
open diamond. The three estimates of Bouwens et al. (2003a) are 
shown by crossed circles and indicate three different complet enes s 
corrections for one sample of objects (see discussion in Section l5.5l . 
The lower panel indicates the star formation density corrected for 
extinction using the Calzetti (1997,2000) reddening law and the 
typical E(B — V) = 0.15 value found by Steidel et al. for their 
z ~ 3 sample and applied by them at z ~ 4 and by Iwata et al. 
at z ~ 5. Note that this assumption of no evolution in the dust 
characteristics of Lyman break galaxies between z ~ 3 and z ~ 6 is 
made primarily for comparison with the work of previous authors. 



is itself evolving lockstep with redshift in the observed 
population (i.e., galaxies would have to be more heav- 
ily obscured at early epochs to explain our results if the 
global star formation rate remains roughly constant). 

Vijh et al. (2003) made a study of the dust character- 
istics of Lyman-break galaxies using the PEGASE syn- 
thetic galaxy templates. They report an ultraviolet flux 
attenuation factor due to dust of 5.7 — 18 at 2 < z < 4 
(c.f., Steidel et al. 1999, who derive a UV attenuation 
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factor of 4.7 using the reddening formalism of Calzetti 
et al. 1994, 2000). Vijh et al. find no evolution in dust 
characteristics of Lyman break galaxies in this redshift 
range. Hence our assumption of no evolution in the dust 
properties of Lyman break galaxies from z « 3 to 2 ~ 6 
(when the Universe was only ~0.55 Gyrs old) may not 
be unreasonable. The lower panel of figure shows the 
extinction-corrected volume-averaged star formation his- 
tory. 

5.3. The Evolution of Lyman-a Emission? 

For the total sample (GOOD-S + GOODS-N), we con- 
firm very few i'-drops as being high-redshift galaxies in 
Lyman-a emission (2 definite and one possible out of a 
total of 10 i'-drop objects observed spectroscopically). 
However, we note that of the 10 z'-drops targetted with 
Keck spectroscopy, two have stellar spectra (GOODS-N 
#1 & #2) and a further two have near-infrared colours 
of low-mass stars (GOODS-N #10; GOODS-S #5), all 
four being unresolved point sources in the HST/ ACS im- 
ages. So, excluding the probable stars, we have strong 
Lyman-a emission in two out of six galaxies, one third of 
our spectroscopic sample (SBM03#1 W T J S ^ = 30 A, this 
paper ; SBM03#3 = 20 A, paper 2). This is in fact 

consistent with the fraction of strong Lyman-a emitters 
at 2 ~ 3: Shapley et al. (2003) report that 25 per cent 
of Lyman-break galaxies have Lyman-a in emission with 
rest frame equivalent widths > 20 A. 

Now we consider whether we would be sensitive to 
emission lines of lower equivalent width in our spec- 
troscopy. Our 5 a sensitivity to line emission is 3 x 
1Q-18 er g S cm -2 g -i ^ or a i_ arcsec extraction and a 

300 km s -1 line- width). Therefore, at the magnitude 
cut-off of our sample {z' AB — 26.5), this corresponds to 

W^rest* > 5A at z « 6. The fact that we do not observe 
any other Lyman-a emission lines, despite our sensitive 
equivalent width limit, is marginally at odds with the 
global properties of the Lyman-break sample at 2 ~ 3 
reported in Shapley et al. (2003): 20 per cent of the Shap- 
ley sample have Lyman-a in emission with < 20 A 
- so we might expect to have seen another emission line 
in our sample at lower equivalent width than the two 
emission lines found. 

There is some evidence for a decline in the fraction 
of Lyman-a emitters at z ~ 6 compared with z ~ 3, 
although the statistics are marginal. If this proves to 
be the case, it could imply that the contamination of 
our total sample by low- 2 galaxies is significant, or per- 
haps that Lyman-a does not routinely escape from z ~ 6 
galaxies. Lyman-a is seen in emission in half of star 
forming galaxies at z ~ 3 (see Steidel et al. 1999). If, 
at 2 ~ 6, we are at the end of reionization then neutral 
gas in the Universe may resonantly scatter this emission 
(the Gunn-Peterson effect) and so lead to a fall in the 
observed fraction of Lyman-a emitters. This effect may 
have been detected by Becker et al. (2001) at z = 6.3 and 
could indicate either that reionization ended later than 
estimated by the WMAP results (Kogut et al. 2003) or 
that pockets of neutral gas remained well after the end 
of the bulk of reionization. 

Maier et al. (2003) describe the Calar Alto Deep Imag- 
ing Survey (CADIS) search for high redshift Lyman-a 
emitters. They also find that bright Lyman-a emitters 



are rarer at z > 5 than suggested by a non-evolving pop- 
ulation, indicating a possible fall in star formation rate. 

This conclusion is supported by recent results from the 
8m Subaru telescope which has now been used for a num- 
ber of searches for high redshift Lyman-a emitters us- 
ing narrow band filter selection. Taniguchi et al. (2003) 
provide a review of recent observations of high redshift 
Lyman-a emitters. They reproduce the Lyman-a emit- 
ter 'Madau' plot from Kodaira et al. (2003) and report 
a fall by a factor of 4 in the global star formation rate 
between z ~ 3 and z ~ 6.6, mirroring the decline we 
from our Lyman-break-selected i'-drop sample. 

5.4. Impact of Evolution of the Luminosity Function or 
IMF on Determining the Global Star Formation 
Rate 

It is fundamental to any discussion of the evolution of 
basic properties such as the comoving volume-averaged 
star formation density that like must be compared to 
like. Since our comparison with the work of previous 
authors is dependent on assumptions of non-evolution 
in both the stellar initial mass function and luminosity 
function of Lyman break galaxies we now consider the 
possible influence on our result of evolution in either of 
these functions. 

Since the UV luminosity-SFR conversion factor extrap- 
olates the total star formation rate given the observed 
flux from blue, high mass stars, it is sensitive to the shape 
of the IMF assumed. In common with previous authors 
we use a Salpeter IMF, however, if the Scalo (1986) IMF 
is used, the inferred star formation rates are a factor of 
« 2.5 higher for a similar mass range. 

Recent work by several authors, (e.g Abel et al. 2002, 
Nakamura & Umemura 2002) has suggested that the ear- 
liest population of stars in the universe (population III) 
has an initial mass function (IMF) that is biased, with 
respect to the Salpeter IMF used here, towards massive 
stars. This implies that if such early stellar populations 
are dominant at z ~ 6 then we are likely to be overes- 
timating the star formation rate (although this clearly 
depends on the detailed shape of the IMF). 

Recent results (Kogut et al. 2003), however, favor an 
end to reionization as early as z ~ 11 (depending on 
the reionization history assumed). By z ~ 6, therefore, 
the era of population III stars may well be at an end. 
The semi-analytic models of Somerville and Livio (2003) 
also suggest that the ionizing flux due to population II 
stars (assumed to have a Salpeter IMF) outstrips that of 
the more massive population III by a factor of ~ 100 by 
2^6. As a result, we consider that use of a standard 
IMF is more appropriate for our calculation at 2^6. 

As described in section |21 our magnitude limit corre- 
sponds to m* for a Lyman-break galaxy at 2 ~ 6 if no 
evolution is assumed in the Lyman-break galaxy lumi- 
nosity function found at 2 ~ 3 and applied at 2 ~ 4 
by Steidel and coworkers (1999). In figure we chose 
to recalculate the reported global star formation den- 
sities of other authors to allow a comparison with our 
SFR> 15 hjQ Mq yr _1 limit and observe a significant fall 
in star formation density when making this assumption 
(although an assumed luminosity fucntion is necessary 
in order to calculate the effective volume of the survey). 
This apparent fall may also arise if the luminosity func- 
tion evolves significantly in this redshift range (a period 
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Fig. 10. — Evolution of the Luminosity Function. In the event 
that there is no fall in the comoving volume averaged star forma- 
tion density between z ~ 3 and z ~ 6 our apparent fall could be 
explained by a variation in the parameter of the galaxy luminos- 
ity function as shown here (crosses). Also shown are the case of a 
moderate fall in the star formation density (asterisks) and for a fall 
in the luminosity function as shown in ngure|5](diamonds). Dotted 
lines indicate the parameter values found at z ~ 3 by Steidel et al. 
(1999). 



between 0.9 and 1.5 Gyr after the Big Bang). 

The galaxy luminosity function ^(L) is conventionally 
parameterized by a Schecter (1976) function of the form: 

*(^ = *'exp(zi) 

Steidel et al. (1999) found no conclusive evidence for evo- 
lution of the rest-frame ultraviolet luminosity function 
between z ~ 3 and z ~ 4 deriving a a faint-end slope 
a = —1.6, a typical magnitude m* = 24.48 and a nor- 
malization of 4>* — 3.5 x 10~ 2 , hj Mpc ~ 3 (in our cos- 
mology) at z ~ 3. If we assume no evolution in a or L* 
then our result corresponds to a 0* =6 = Q.l(f>*_ 3 . 

Figurc^Hillustrates the locus of parameter values that 
would produce our results if we instead assume that the 
comoving star formation density remains unchanged be- 
tween z ~ 3 and z ~ 6 . Also shown are locii correspond- 
ing to a moderate fall in the star formation density and 
for a fall as shown in figure 

As can be seen, if the star formation density is in fact 
unchanged from that at z ~ 3, then our results may be 
explained either by a fall in the characteristic luminosity 
of Lyman break galaxies at z ~ 6 together with a corre- 
sponding rise in object number density (L* =6 = 0.5L* =3 , 
a = —1.6, 4>l =e = 2</>* =3 ) or by a steepening of the faint 
end slope (a = -2.0, L* =6 =L* =3 , </>* =6 = 0* =3 ) or by 
some combination of the two. In practise, either extreme 
falls outside the expected range of values for these pa- 
rameters at z ~ 6. 

5.5. Comparison with Previous Studies 

In recent months a number of groups have made studies 
of Lyman break galaxies and star formation histories at 
z > 5 utilizing deep space or ground based imaging. In 



this section we compare these studies and contrast their 
results with those presented above. 

Multi-wavelength ground-based imaging and the pho- 
tometric redshifts that may be derived from it are valu- 
able tools in the identification of high redshift objects 
but are limited by the resolution and magnitude limits 
that can be reached. Fontana et al. (2002) presented a 
sample of z > 4.5 candidates selected at z' AB < 25.0 from 
imaging with the Very Large Telescope (VLT). They find 
13 high redshift candidates including 4 in their highest 
redshift bin of 5 < z < 6 in a total area of 29.9 arcmin 2 
and observe an order of magnitude decrease in the UV 
density between z ~ 4.5 and z ~ 6. However, correct- 
ing for their selection function and absolute magnitude 
limit in each redshift bin they report no decrease in the 
comoving global star formation rate. This is a relatively 
small survey and the authors note significant variance 
between their two survey fields. In figure^] we show star 
formation densities derived using our prescription from 
the UV luminosity densities reported by Fontana et al. 
for their magnitude-limited 'minimal' sample of objects 
which they are confident lie at high redshift. As can be 
seen their results are in good agreement with our own 
although based on very small numbers of objects. 

HST imaging such as that studied in this paper allows 
such surveys to be pushed to fainter and more distant ob- 
jects. At z > 4.5 Lehnert and Bremer (2003) presented 
a sample of objects selected on the basis of their large 
(R — 7)ab color. Their number counts at Tab < 26.2 
are slightly lower than expected if no assumption in the 
luminosity function is assumed but not significantly so 
given the effects of completeness corrections and cosmic 
variance. Similarly Iwata and coworkers (2003) consider 
a catalog of candidate z ~ 6 objects selected primarily 
on the basis of large V — I colors. The star formation 
density estimated by Iwata et al. from their survey in- 
dicates a slight fall in global star formation density in 
the range 4 < z < 5 although their result could also be 
interpreted as consistent with a constant star formation 
density at z > 2. Capak et al. (2003) comment on the 
selection criteria in both these papers and suggest that 
both groups may have underestimated the contamination 
in their sample due to interloper galaxies at lower red- 
shifts. If this is the case, then the star formation density 
derived by Iwata et al. for z ~ 5 (shown in figure |5J) may 
in fact represent an overestimate. This would support 
our observations of a fall in global star formation density 
with redshift. 

Yan et al. (2002) also used deep HST imaging to iden- 
tify 30 z ~ 6 candidates in an area of 10 armin 2 to a 
depth of z' AB = 28.3 and with a median magnitude of 
z' AB = 21 A. They estimate their possible contamina- 
tion from cool dwarf stars and elliptical galaxies as 7 ob- 
jects (23 per cent), a lower fraction than that estimated 
by most other authors. None of the objects identified 
by these authors would be selected at our limiting mag- 
nitude, consistent with our expected 0.2 ± 0.1 objects 
in a 10 arcmin 2 field. In common with Bouwens et al. 
(2003a), we note the difficulty of identifying objects so 
faint with certainty given the total exposure times quoted 
by Yan et al. Nonetheless, we note with interest the sug- 
gestion made by these authors that their number counts 
at faint magnitude suggest that the faint end slope, a, 
of the galaxy luminosity may be as steep as a = — 2 (see 
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section IQ1 for further discussion of this). 

Finally a selection of z ~ 6 candidate objects made 
from two deep fields observed with the HST /ACS in- 
strument has been presented by Bouwens et al. (2003a) 
They find a total of 23 objects to a magnitude limit of 
z' AB = 27.3(6er), somewhat fainter than our limit, and 
use their data (selected with a similar (i' — z')ab > 1.5 
cut to the sample in this paper) to make three estimates 
of the star formation density at this redshift. The third 
estimate, presented by Bouwens and coworkers as a strict 
lower limit, is largely consistent with the results pre- 
sented in this paper (given the different survey depth) 
and is derived using a similar methodology and assump- 
tion of no evolution in the intrinsic characteristics of Ly- 
man break galaxies between z ~ 3 and z ~ 6. However, 
it is notable that their sample includes only two objects 
that would be selected in this paper, broadly consistent 
with the 1.1 ± 0.5 objects we would expect to observe in 
their 46 arcmin 2 survey area. Their sample also includes 
one very luminous object which substantially increases 
their final result. In addition Bouwens et al. present 
two incompleteness-corrected estimates which make use 
of the 'cloning' procedure described in Bouwens et al. 
(2003b). This procedure projects a sample of V-drop se- 
lected objects (selected to lie at z ~ 5) to z ~ 6 taking 
into account the effects of k-corrections, surface bright- 
ness dimming and PSF variations. This sample is then 
used to assess their survey incompleteness and selection 
function, leading to star formation densities some 2.7- 
3.8 times higher than their third estimate. This proce- 
dure takes into account the natural variability in size, 
shape and surface brightness of galaxies but in doing so 
makes the implicit assumption that galaxies at z ~ 6 
have similar distributions in these properties to galax- 
ies at z ~ 5. In addition, their first (higher) estimate 
is based on comparison with their earlier work at lower 
redshifts (Bouwens et al. 2003b) which is based on a rel- 
atively small observation area which substantially influ- 
ences their final result. 

The observed trend in the global star formation density 
at z > 3 is still uncertain and the evidence mixed. Al- 
though this work suggests a fall in star formation density 
based on the bright end of the luminosity function, other, 
deeper, surveys suggest that it is possible that there is 
no fall in the integrated star formation rate. This sug- 
gests that the shape of the luminosity function may well 
be evolving over this redshift interval, as discussed in 
section 15.41 The results of Bouwens and coworkers in 
particular, however, illustrate the complexities involved 
in fully assessing the completeness and selection func- 
tion of a photometrically derived galaxy sample at z ~ 6 
and the necessity of making any completeness corrections 
assumed completely transparent in order to allow com- 
parison of work in this field. 

5.6. Comparison with the results of the GOODS team 

During the refereeing process of this manuscript a num- 
ber of papers written by the GOODS team (described in 
Giavalisco et al. 2003a), and based on three co-added 
epochs of GOODS data, have become public. 

A colour selection similar to that in this paper and that 
of Bouwens et al. (2003) was performed by Dickinson et 
al. (2003) in order to identify z w 6 candidates with a 
slightly bluer colour cut criteria of (i' — z 1 ) ab > 1-3 which 



may leave them susceptible to increased contamination 
from lower redshift galaxies. Rather than apply a strict 
magnitude limit, these authors have chosen to apply a cut 
based on object signal-to-noise in their co-added images. 
They have also limited their survey area to that observed 
in all three epochs (approximately 80% of the total area) 
and limited their sample to resolved objects. As a re- 
sult, a large fraction of the z ~ 6 candidates described 
in this paper and in paper 1 do not form part of their 
survey. Dickinson et al. identify only five candidate ob- 
jects with signal-to-noise greater than 10 in the combined 
GOODS-N and GOODS-S survey of which three objects 
(in the GOODS-S field) were also selected in paper 1 and 
the remaining GOODS-S object falls outside our colour 
and magnitude cuts. The remaining candidate, NiDOOl, 
has colours which appear to vary significantly between 
epochs and averages across 4 epochs as too faint to ap- 
pear in our selection. These authors also present a deeper 
sample of galaxies with signal to noise > 5 although they 
acknowledge that this is likely to contain significant con- 
tamination. They note that the observed number density 
of bright galaxies is significantly lower than that expected 
from the lower redshift luminosity functions for Lyman 
break galaxies, supporting our conclusion that the lumi- 
nosity function of these objects must evolve (either in 
shape or normalisation) between z ~ 3 and z ~ 6. 

Giavalisco et al. (2003b) utilise the Dickinson et al. can- 
didate selection at z ~ 6, together with similar colour 
selections at lower redshift, to comment on the evolu- 
tion of the rest-frame UV luminosity density (and hence 
star formation rate) with redshift. These authors have 
attempted to correct the low signal-to-noise Dickinson 
et al. sample for contamination and incompleteness and 
work to a faint magnitude limit. This yields a star for- 
mation history that declines only slowly with redshift if 
evolution in the luminosity function is modelled using a 
X 2 approach (which assumes a priori a = —1.6) and a fall 
if the V e f f method using the z = 3 luminosity function 
is adopted. While the x 2 approach has the advantage 
of simultaneously fitting the luminosity function and al- 
lowing an estimate of the sample volume, and dows not 
reuire an assumed lumosity function, it is not clear that 
the number counts of z ~ 6 candidates are sufficient to 
support this approach (Dickinson et al. only report five 
objects with confidence). However, as Giavalisco et al. 
comment, their result is largely in agreement with that 
presented here and in paper 1 , given the observed decline 
in number density of bright galaxies at z ~ 6. Their re- 
sult of an almost constant luminosity density at z > 3.5 
is highly dependent on low-significance faint end counts 
and on the corrections applied to these. It seems possible 
that the number counts of faint objects will indeed show 
evolution in the shape of the luminosity function between 
z = 3 and z — 6, yielding a flat star formation history 
at high redshift but the data are not currently available 
to determine this. As Giavalisco et al. comment, the 
luminosity density measurement at z ~ 6 clearly needs 
to be revisited with deeper data, which may be possi- 
ble with the forthcoming release of the HST Ultra Deep 
Field (UDF). 

6. CONCLUSIONS 

In this paper we have presented a selection of objects in 
the Hubble Deep Field North with extreme i' — z' colors 
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and identify a number of good z ~ 6 candidates. From 
these we have calculated an estimate of the global star 
formation density at this redshift. Our main conclusions 
can be summarized as follows: 

1. The surface density of i'-drops with [i 1 — z')ab > 
1.5 and z' AB < 25.6 in the GOODS-N is 0.03 
objects/arcmin 2 , similar to that in GOODS-S re- 
ported in paper 1. This indicates that cosmic vari- 
ance is not the dominant uncertainty in measur- 
ing the star formation density at z ~ 6 from the 
GOODS EST /ACS images although it is possible 
we may have the first suggestion of a group at 
z = 5.8 in the GOODS-S field. 

2. We measure the star formation rate from the rest- 
frame 1500 A flux in z ~ 6 galaxies with ob- 
served star formation rates > 15 Mq yr _1 (un- 
corrected for extinction). We detect a decline in 
the comoving volume-averaged star formation rate 
between z ~ 3 to z ~ 6 for resolved objects at this 
bright limit, consistent with the results reported in 
paper 1. The observed fall if all i-drops (includ- 
ing unresolved objects) in the field are considered 
is less pronounced but this is likely to signifcantly 
overestimate the star formation density at z ~ 6. 

3. We have considered the possible effect of evolution 
in the shape and normalization of the luminosity 
function; if the total comoving star formation den- 
sity (including galaxies well below our flux limit) is 
in fact the same as at z ~ 3, then the luminosity 
function must change significantly (from z = 3 to 
z = 6, L* drops by a factor of 2 if <jf and a are 
constant, or alpha changes from —1.6 to —2.0 if cf>* 
and L* remain constant at z > 3). 

4. Ultra-deep Keck/DEIMOS spectroscopy has been 
obtained for half our sample. The long-slit spec- 
trum of GOODS-N object 4 shows a single emission 
line at 8807 A, offset by 2" , although it is possi- 
ble that this is associated with a nearby galaxy at 
lower redshift (not an i'-drop). If it is spatially- 
extended Ly-a emission then the redshift would be 
z = 6.24, although a lower-redshift interpretation 
is possible. We detect continua from two sources, 
with low signal-to-noise, but these may be consis- 
tent with late-type stars. Two fainter sources have 
no detectable flux in the spectroscopy. 

5. For the total sample (GOODS N+S), we confirm 
very few i'-drop galaxies in Lyman-a emission (2 



or 3 out of ~ 15). This could be that the contami- 
nation of low-z galaxies and stars is significant, or 
perhaps that Lyman-a does not routinely escape 
from z ~ 6 galaxies. 



6. We have analyzed surface brightness selection ef- 
fects: although undoubtedly these play a role, in 
fact the vast majority of our candidates are very 
compact (as would be expected for star-forming re- 
gions), so this seems not to be a significant effect. 
We correct our survey volume for the dimming in 
apparent z'-band magnitude with redshift by using 
a V/Vmax "effective volume" approach, introducing 
a correction of about a factor of two. 
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